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■ Using the symmetry properties of the three-anyon spectrum, we obtain ex- 

actly the multiplicities of states with given energy and angular momentum. The 
[ results are shown to be in agreement with the proper quantum mechanical and 

semiclassical considerations, and the unexplained points are indicated. 

> 
• 1— 1 

X 



email: gezin@gluk.apc.org 



It is well known that the quantum mechanical spectrum of N non-interacting bosons 
or fermions can be obtained given just the single-particle spectrum, but for anyons 
this does not hold because the iV-anyon problem is essentially many-particle. For 
arbitrary iV there are two classes of exact solutions [[1], |2], |3|, Q but their relative number 
decreases rapidly with TV increasing. The only case in which one can proceed with the 
exact analysis more or less far ahead is that of N = 3. In the previous work || we 
have shown that it is possible to calculate exactly all the degeneracies in the three- 
anyon spectrum using certain symmetry properties of the latter. Here we will carry 
out an analogous calculation, taking into account in addition the angular momentum. 
Our present consideration will allow us to shed at least some light on the problem of 
quantum mechanical description of anyonic spectra, which at the moment is far from 
being closed. 

As it was done earlier, we consider the problem of three non-interacting anyons 
in a harmonic potential, with the particle mass and the frequency set to unity. The 
Hamiltonian is H = Y^=±Hj with the one-particle Hamiltonian Hj = |(— Aj + r 2 ). 
The single-particle state is uniquely determined by the two quantum numbers - energy 
E which may equal 1,2,..., and angular momentum L = —(E — 1),—(E — 3),...,E — 
3, E — 1. Formally, the number of single-particle states with energy E and momentum 
L is given by 

f r(E-L,2) if \L\ < E-l 
9i(E,L) = l (1) 
[ otherwise 

where r(a, b) is the remainder of the division of a by b. (Here and further E is implicit 
to be a positive integer and L an integer.) In what follows it will be also convenient to 
use, along with E and L, the numbers i + = (E + L — l)/2 and t~ = {E — L — l)/2. 
Obviously, any pair of non-negative integers (i + ,i~) corresponds to exactly one state. 

To calculate the multiplicities g 3 B(E, L) and g 3 F(E, L) in the relative motion spec- 
trum of three bosons and three fermions, respectively, we start from the multiplicities 
gs(E, L) corresponding to Boltzmann statistics, then come to g 3 B(E, L) and g3F{E, L) 
for the full spectrum, and afterwards separate the center-of-mass motion. 



The calculation is simplified by using 



(2) 



since C^ 1 = if + if + if with the Ef's as defined earlier, the Boltzmann degeneracy is 
gs(E, L) = S3(C + )Ss(C^) where Ss(N) is the number of ordered triples of integers the 
sum of which is N. One has S 3 (N) = ^(N + l)(N + 2) and consequently 



g 3 (E, L) = 1 [(E 2 - l) 2 - 2(E 2 + 1)L 2 + L A 



(3) 
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Now, in Boltzmann count one bosonic state is taken six times if all the particles occupy 
different single-particle states, three times if two are in the same state but the third is 
in another one and one time if all are in one and the same state. Therefore 



g 3B (E, L) = ± [g 3 (E, L) + 3d 3 (E, L) + 2t 3 (E, L)] 



(4) 



where d 3 (E, L) is the number of states with two particles in the same state (the third 
one may be or not be in that state) and t 3 (E, L) is the number of states with all three 
particles in the same state. For fermions, correspondingly, 



g 3F (E, L) = - [g 3 (E, L) - 3d 3 (E, L) + 2t 3 (E, L)\ . 
o 



(5) 



By virtue of the aforesaid, 

d 3 (E,L)=Q 3 (£+)Q 3 (£-) (6) 
with Q 3 (N) the number of ordered pairs (£i,£ 3 ) such that 2£± + £ 3 = N {£2 = £i), 



N 
~2 



+ l = -[JV-r(iV,2)]+l, 



Qs(N) 

where [n] stands for the entire part of n. Finally, 

t 3 (E,L) = d{E,3)d(L, 3), 
where we have introduced the "multiplicity function" 

1 if r(a, b) = 

otherwise 
Explicit expressions for g 3 B and g 3 F read 



(7) 



(8) 



d(a, b) 



(9) 



93b(E,L) 



384 



(E 2 - L 2 ) 2 + WE 2 - UL 2 + 2AE + 13 



+ ^-[2M + M- + M + (L -E-1)-M-(L + E + 1)} + \d{E, 3)d(L, 3), 
lb 3 



(10) 



93f(E,L) 



1 

384 



(E 2 -L 



2\2 



UE 2 + 10L 2 - 2AE - 11 



— [2M + M- + M + {L — E — 1) — M-{L + E + 1)1 + -d(E, 3)d(L, 3), 
16 3 



where M ± = r(£ ± ,2) , C^ 1 are defined by (§), and E and L should be of different 
parity. 
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In order to separate the center-of-mass motion, note that 



g 3S (E,L) = J2 E hs(E-e,L-l) (12) 

e=l l=-e+l 

for S = B, F ( e and I are the center-of-mass energy and angular momentum, respec- 
tively). Eq.(|l2"D implies 

g 3S (E, L) = g 3S (E + l,L) + g 3S (E- 1, L) - g 3S (E, L - 1) - g 3S {E, L + 1). (13) 

Substituting (|l(J) and (|Tl~D , one gets 

g 3{B /F)(E, L) = ^(E 2 - L 2 ) + i + 1, 3)d(L, 3) + d(E - 1, 3)d(L, 3) 

-d(E, 3)d(L -1,3)- d(E, 3)d(L + 1, 3)] ± \r(E - L, A)r(E + L, 4), (14) 

8 

the upper/lower sign referring to B/F ; here, E and L have to be of the same parity. 
In the table below, the values of g 3 s{E, L) and g 3 p{E, L) for all E < 12 are listed. 



Now, our main goal is to investigate how the spectrum interpolates between the 
bosonic and fermionic ones. Anyonic wave functions satisfy the interchange conditions 
Vjk^ = exp(z7T<5)\l/ for each pair j, k , where Vjk is the operator of anticlockwise 
interchange of particles j and k and S is the statistical parameter; there is a continuous 
transition from bosons to fermions as S goes from to 1 . Having turned an iV-anyon 
system by one complete revolution in an anticlockwise direction, one will have the 
wave function multiplied by exp[z7riV(iV — 1)5] . Consequently, the possible values 
of angular momentum for anyons are L = jV ( j ^" 1 ) £-|-integer. In our case N = 3 one 
gets L = 35+integer . Thus, a state with angular momentum L at Bose statistics 
interpolates to that with angular momentum L + 3 at Fermi statistics. As for energy, it 
is known || [7], || that possible values of the difference E Fermi —E Bose are +3, +1,-1, —3. 
In the first and the last cases the 5 dependence of E is linear and the states may 
be found exactly (in Ref. we called them "good" states), in the other two ones 
("bad" states) this dependence is non-linear and not found exactly (see Ref. for 
an interesting hypothesis concerning the mentioned dependence). To summarize, all 
bosonic states with the quantum numbers (E, L) fall into four classes according to 
the quantum numbers of fermionic states to which they interpolate; those can equal 
(E + 3, L + 3), (E+1,L+ 3), (E - 1, L + 3), or (E - 3, L + 3). (It is implicit of course 
that in the subspace of states with the same {E,L) one chooses the "correct" ones in 
the same way as in perturbation theory with degenerate states; it can always be done 
since E and L themselves are good quantum numbers for anyons.) 

We will use r n (E, L) to denote the number of states which come from (E, L) at Bose 
statistics to (E + n, L + 3) at Fermi statistics. (Since it has been shown by numerical 
calculations 0, that the energies of different "bad" states interpolating between the 
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same bosonic and fermionic ones are in general different for fractional 5 , one has for 
anyons, strictly speaking, "numbers" of such states rather than "multiplicities".) As 
n can take four values, one needs four equations for each (E, L) to determine those 
numbers. The "law of conservation of states" at bosonic and fermionic points reads 

r +3 (E, L) + r +1 (E, L) + f~\E, L) + r~ 3 (E, L) = ~g 3B (E, L), (15) 

r +3 (E — 3, L — 3) + f +1 (E — 1, L — 3) 
+ r\E + 1, L - 3) + r~ 3 (E + 3, L - 3) = g 3F (E, L), (16) 
respectively. Two more equations follow from the known symmetry properties of the 



spectrum. First, in perturbation theory it is easy to establish [|10| , |TT|| that at Fermi 
statistics, two states with opposite angular momenta have opposite values of slopes, 
that is, derivatives {dE/d5)s=\- Since the linear behavior of the "good" states is exact, 
this implies the equality 

f +3 (£-3,L-3) =f- 3 (E + 3,-L-3) (17) 

(see Fig.l). Second, there is the supersymmetry property. Namely, it was pointed out 
by Sen |T2[ that there exists an operator Q which annihilates some of the "good" states 
but acting on a "bad" state with statistical parameter 5 always produces another "bad" 
state with same energy and with statistical parameter 1 + 5. It is straightforward to 
show that [L, Q] = 2Q so that a state coming from (E, L) at 5 = to (E + 1, L + 3) at 
5=1 turns under Q into that coming from (E, L+2) at 8 = 1 to (£'+1, L+5) at 5 = 2. 
Now, parity transformation turns the latter into a state coming from (E + 1, — L — 5) 
at S = to (E, —L — 2) at 5 = 1 (Fig.2). Therefore the last equation is 

f +1 (E,L) = f-\E + 1,-L - 5). (18) 

The four equations (|T5"l) -(|l"8l) and the expression ([14]) for g 3 B and would suffice to 
calculate r n (E, L). It seems, however, hardly possible to obtain closed-form expressions 
for them in a straightforward manner as it is not clear in which form they should be 
searched for. Instead, it is easy to solve ([0f)-(ll8D numerically for low- lying levels and 
it turns out to be possible to pick up the regularity in the numbers. 

The results for r n (E, L) for E < 12 are summarized in Tab. 2. The general tendency 
is clear and can be said to coincide with what one could more or less expect, but to 
move further it is useful to involve the concept of towers |], [TTJ, jT2J. Again following 



Ref. [12 1, we recall that there exists an operator K_ such that [H, KJ\ = —2K_ and 
[L, K-] =0 — therefore it either annihilates a common eigenstate of H and L or lowers 
its energy by two units without changing its angular momentum. This means in turn 
that all the states fall into "towers" descending along which is realized by K_ until 
one reaches a bottom state for which = 0. Therefore any (E, L) state either is a 

bottom state or has its correspondent {E — 2, L) state obtained from it by action of 
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. Denoting by b n (E, L) the number of bottom states with the same characteristics 
as in r n (E, L) , we have 



b n {E, L) = f n (E, L) - r n {E - 2, L), 



(19) 



and correspondingly 



[E/2] 

r n (E,L) = J2 b n (E-2k,L). 

k=0 



(20) 



Thus, instead of counting (and finding) all states it is sufficient to count (and find) 
only bottom states. The values of b +3 (E,L) and b +l (E,L) are displayed in Tabs. 3-4. 

Here at last the regularities are obvious. The formal expressions read 



b +3 (E,L) 



E + V 




~E + L~ 




~L~ 


[ 6 . 




[ 4 


+ 


.2. 




- 2 






for — 




< L < E 


-2; 




3 







l-d(E,6) 



b +1 (E,L) 



~ E*~ 




[ 6 . 


+ 



L-L* 1 - r(E,2) 



4 



+ 



r(E*,6) 



[l-2r(E,2)} 



d(E* -2,6)[r{L- L*,4)-l] 



for < L < E - 6, 



where 



E* = E-3r{E,2), L* = 2 



E* 
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2 + r(E,2). 



From these, one obtains the formulas for f n (E, L) immediately by applying 

f L 2 + QL + 5 1 - r(L, 2) d(L, 3) 



r +3 (E,L) = { 



12 



4 



for < L < 2=1, 



—E 2 + 6EL - 5L 2 + 12E - 12L d(L, 3) - d(E, 3) 



+ [1 -2r{E,2 



48 

d(E-L -2,4) 



4 



for < L < E — 2; 



(21) 



(22) 
(23) 



(24) 
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' E 2 + QEL + 9L 2 + 16E + 48L + 48 r(E - L,A) d(E- 1,3) 



r + \E,L) 



48 



+ 



for =f=2 < L < —1, 



E 2 + — 15L 2 + 16E — 72L — 96 r(E - L,A) d(E - 1,3) 



48 



+ 
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+ 



for — 1 < L < 



3 ' 



E 2 -2EL + L 2 -AE + AL r(E — L, A) 



16 



+ 



8 



for < L < E — 6. 



(25) 



Formulas (|2T|) , ([Zzj) , (|24|) , (fig) are exact and are the main result of this paper. In 
these formulas, E and L should be of the same parity, i.e. r(E — L, 2) = ; otherwise, as 
well as if L does not fall in any of the ranges specified in the formulas, the corresponding 
multiplicities vanish. Formulas for 6 _1 ,f _1 and 6~ 3 ,f -3 follow immediately from the 
obtained ones upon applying (|17D and (|T^). 



Finally, summation over L yields the total numbers of states with given energy and 
slope: 



(E) = 2l6 + 9E " + [42 " 2) " 24rf(E ' 3)] E + [ /r(£ *' 6) " 81r(jB ' 2) 1 } 



(26) 



r +1 (E) = — {2E 3 + 3E 2 + [2Ad(E -1,3)- 18] E + [Ar(E*, 6) - 21r{E, 2)]} , (27) 

where fo — 0, ji — 88, f\ = 56, and E* has been defined by (p3|); these are exactly the 
expressions obtained in Ref.[Q (where slightly different notations were used). 

Let us now discuss the obtained results. For clarity, the plots of f n (E, L) for 
E = 100 are displayed in Fig. 3. We have learned that for a given E , the states of each 
of the four classes exist for not all values of L allowed for that E . Tab. 5 shows the 
lowest and the highest values of L for which the corresponding states exist, as well as 
the values for which the numbers of such states are maximal. 



Since all the states are solutions of the Schrodinger equation, all these values should 
in principle be deduced directly from this equation. Producing them in this way is an 
interesting and apparently difficult problem, and perhaps at least some points of its 
solution, if available, could be used also to understand the structure of the iV-anyon 
spectrum. What we can do at the moment is to explain why for L = E — 2 and 
E — A , as the table shows, there exist only the (+3) states. In the context of our 
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previous considerations, this follows immediately from (|ITD and ( |i~8|) by noticing that 
f +l (E, E — k) = f^ 1 {E + 1, —E + k — 5) does not vanish only for k > 6, because there 
are no (E, L) states with L < —E + 2 , and analogously for f -1 and f~ 3 . However, 



a purely quantum mechanical proof of this fact can be given 0, 13], which we will 
describe here. 



Choosing as independent coordinates the one of the center of mass Z = [z\ + z 2 + 
Zz)/V% {zj = Xj + iyj is the complex coordinate of j-th particle) and two relative 
z 2 )/\/2 , Z23 = (z2 — z%)j\pl and searching for the wave function in 



ones Z\i 

the form ^ = ^ cm X ex P 

H = H(t) + H(2) where = z 12 d 12 + z* l2 d* l2 + -223^23 + zh^h + 2 , H, 



— -Az\Z\ + Z2Z\ + Z3Z3 



one has for the relative Hamiltonian 



(2) = -2(012% + 

<9 2 3<9 23 ) > fyk = d/dzjk ( H acts on x), and for the relative angular momentum 
L = Z\ 2 d\2 — z \2®\2 + ^23023 — Z 23^23- The function x is to be searched for as a linear 
combination of functions of the form 



ho ho : Z93 I 



12 '12 1 '23 '23 i '31 '31 



/si h 



Zl2) h2 [- 



12 / 12 (z 23 ) i23 U 



(28) 



^31 = —Zyi — Z2z) where {•••}+ means symmetrization over Zj's and the equalities lj 



ik 



Ijk — 5+integer should fulfil. The function (EHf) then satisfies the anyonic interchange 



conditions VjkX = ex P(^ 7r ^)x> an d the coefficients in the linear combination should be 
chosen so that x be non-singular and satisfy the equation H\ = E\. Obviously, 



L 



ll2 ll2 ', hs h3 ! ^31 ^31 



(jl2 — ^12 + ^23 ~~ ^23 + ^31 — ^31 



^12 ^12 ! ^23 ^23 j ^31 ^31 



H, 



(i) 



^12 ^12 ! ^23 ^23 ! ^31 ^31 



(/l2 + I 



12 T I23 + ^23 + ^31 + ^31 + 2 



^12 ^12 ! ^23 ^23 ! ^31 ^31 



and H, 



(i> 



1 12 ^12 ; ^23 ^23 ; ^31 ^31 / consists of pieces of the form 



(29) 



(30) 



• Ijk 1 • • • Imn 1 



The states which at 5 = have the maximal angular momentum L = E — 2 are 
those with all Ijk vanishing, i.e. of the form |Z 0;m0;n 0) + (with integral l,m,n); the 
corresponding anyonic states \l + 5 0; m + 5 0; n + 5 0) always are stationary because 
H(2) annihilates them, and their energy is E = Esose + 35 (where Eboss = I + m + 
n + 2). Further, to construct a state which has L = E — 4 at 5 = one should start 
from \l + 6 1; m + 6 0; n + 8 0) + = Acting on this with H^) yields a sum of the 
terms of the form \p + 5 0; q + 5 0; r + 5 0) + , which are annihilated by H^) ■ If all 
of them are non-singular, then a stationary state can be built. (One has ifm |1) = 
E\l) , H { 2)\1) = |2) , # ( i)|2) = (E-2)\2) , and H {2) |2) = 0, consequently 
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H 

be 



H) + i|2) 



E 



|1> H" 2 |2> 
5 0), 



Now, the only potentially singular term in |2) could 

with non-negative integral q and r . However, such 
1 0; q 0; r 0) , can easily be shown to be non- 

This completes 



1 + S 0; q + S 0; r - 
function, which equals (212-223-231} 

singular: the singularities cancel out due to symmetrization 
the proof of the fact that all bosonic states with E = L — 4 , as well as those with 
E = L — 2, have their correspondent anyonic states with the slope +3. 
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At present it is unknown how the other features of the obtained results could be de- 
rived from such quantum mechanical considerations. However, the qualitative picture 
can be justified using semiclassical arguments. In the two-anyon problem, there is only 
one relative coordinate Z12 , the change of which is governed by the oscillator equa- 
tion. The corresponding classical trajectories are ellipses, and the anyonic interchange 
conditions demand that the relative angular momentum be L = 5 + 21 , where the 
integer I determines the sign of L at 5 = . Thus, as one starts increasing 5 from zero, 
one should expect linear increase of the energy of those states in which the relative 
vector, in classical terms, rotates anticlockwise (I > 0) or does not rotate (/ = 0) , and 
linear decrease for those in which it rotates clockwise. This is indeed the case: The 
semiclassical description of the two-anyon spectrum turns out to yield the exact values 
of levels |14] , ||. For N anyons, there are N ^~^ relative vectors, and the analogous 
consideration would yield linear dependences with the slopes E(1)—E(0) = N ( N ~V —2s 
with s the number of the relative vectors rotating clockwise. The set of the slopes is 
correct [|T5|], which encourages us to use this picture for a qualitative analysis (although 
we do not argue that it is good in other senses, for example that it gives the correct 
multiplicities). So, for a state to have the slope +3 , neither of the relative vectors 
should rotate clockwise. Certainly it cannot be so if L < . For L = a (+3) 
state can be only realized as a radial excitation of the bosonic ground state; indeed, 
r +3 (E, 0) = 1 (for even E ). Further, if we choose a state with a given L "at random", 
then the more is L , the more vectors "in average" rotate anticlockwise. Therefore the 
states with L close to maximal L = E — 2 have at most the slope +3 , and with L 
decreasing, states with successively decreasing slopes emerge, pass the point of their 
maximal number, and vanish, just as one observes in Fig. 3. Apparently the picture is 
analogous for arbitrary N. 

To summarize, we have determined exactly the multiplicities of states with given 
energies and angular momenta in the spectrum of three non-interacting anyons in a 
harmonic well and shown that the results are in conformity with certain quantum 
mechanical and semiclassical considerations. Generally speaking, these results should 
follow directly from the Schrodinger equation, but the concrete procedure of deriving 
them this way is still to be found. 

I thank G.M.Zinovjev for his attention to my work and constant support, and 
Diptiman Sen for stimulating discussions. This work was supported, in part, by a 
Soros Foundation Grant awarded by the American Physical Society. 
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Tab. 2. The values of r n (E, L). Zeros are not written down for clarity. 
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Tab.3. The values of b +3 {E, L) . 
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Tab.4. The values of b +1 (E, L) . 
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Tab. 5. The lowest and highest values of L for which the states belonging to r n (E, L) 
exist, and the value for which the number of such states is maximal. 
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Figure Captions 

Fig.l. Two "good" states with opposite angular momenta in the fermion limit. 
Fig. 2. Two "bad" states related by (supersymmetry+parity) transformation. 
Fig.3. The plots of f n (E, L) for E = 100 . 
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